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1.  A  non-linear  theory  for  the  sea  response  at  the  bow 
of  ships  pitching  and  heaving  in  head  seas  has  been 
systematically  derived. 

2.  This  theory  has  been  reduced  to  non-linear 
computations  of  the  free  surface  in  the  transverse  plane,  as 
a  function  of  time  (2D&T) . 

3.  The  integral  equations  describing  the  motion  in  the 
transverse  plane  have  been  formulated. 

4.  The  problem  at  the  interaction  of  the  free  surface 
and  the  ship  hull  was  successfully  dealt  with  (this  has  been 
a  difficulty  for  some  other  investigators) . 

5.  A  computer  method  for  the  solution  of  these 
integral  equations  in  time  has  been  formulated,  and  a  code 
prepared. 


6.  Calculations  have  been  successfully  made,  see 
Enclosures  A  &  B. 

7.  The  method  is  able  to  predict  the  history  of  water 
level  along  the  bow,  including  bow  emergence,  slamming,  and 
the  height  of  water  over  the  deckline.  These  are  strongly 
effected  by  non-linear  processes. 

Referencaa 

1.  Annual  Letter  Report  for  FY  90,  from  Prof.  Marshall  P. 
Tulin,  UCSB  to  Dr.  James  Fein,  ONR. 

2.  Annual  Letter  Report  for  FY  91,  from  Prof.  Marshall  P. 
Tulin,  UCSB  to  Dr.  James  Fein,  ONR. 

3.  Annual  Letter  Report  for  FY  92  from  Prof.  Marshall  P. 
Tulin,  UCSB  to  Dr.  James  Fein,  ONR. 

4.  Annual  Letter  Report  for  FY  93  from  Prof.  Marshall  P. 
Tulin,  UCSB  to  Dr.  James  Fein,  ONR. 

5.  OEL  Report  90-53,  "Prediction  of  Deck  Wetness;  A 
Theoretical  Development,"  by  H.  Maruo,  August  1990. 

6.  OEL  Report  91-61,  "Prediction  of  Deck  Wetness  on  Ships 
in  Motion  Through  Waves  at  Finite  Forward  Speed  by  Means 
of  a  Nonlinear  Version  of  the  Slender  Ship  Theory,"  by 
Hajime  Maruo  and  Wusheng  Song,  October  1991. 

7.  OEL  Report  92-76,  "Treatment  of  Corners  on  the  2-D 
Boundary  and  Numerical  Examples  of  Nonlinear  Deck 
Wetness,"  by  Wusheng  Song,  September  1992. 


Final  Report:  N0001490J1245 
Page  3 


8.  OEL  Report  94-104,  "Nonlinear  Analysis  of  Bow  Wave 
Breaking  and  Deck  Wetness,"  by  Hajime  Maruo  and  Wusheng 
Song,  March  1994. 

9.  Song,  W.  and  Maruo,  H..  "Bow  Impact  and  Deck  Wetness: 
Simulations  Based  on  Nonlinear  Slender  Body  Theory, " 
Proceedings  of  the  1993  ISOPE  Conference,  Singapore, 
June  1993. 

10.  Tulin,  M.  and  Wu,  Ming.  "Non-linear  Bow  Waves  on  an 

Inclined  Wedge, "  Proceedings  of  the  1994  Workshop  on 
Waves  and  Bodies,  Kyushu,  Japan. 

11.  Maruo,  H.  and  Song,  W.  "Non-linear  Analysis  of  Bow 

Breaking  and  Deck  Wetness  of  a  High  Speed  Ship  by  the 
Parabolic  Approximation, "  Proceedings  of  the  20th  ONR 
Symposium  on  Naval  Hydrodynamics,  Santa  Barbara,  August 
1994  (to  appear) . 

12.  Wu,  Ming.  Doctoral  Dissertation  (to  appear  in  early 
1995)  . 


5 


Accesion  For 


NTIS  CRA&I 
DTIC  TAB 
Unannouiiced 
Justification 


□ 

□ 


By, . . 

Distribution  / 


Availability  Codes 


Oist 


Avail  and/or 
Special 


Knclosore  X 


Summary  of  Material  Prepared  at  UCSB  During  1992-93 

By 

Ming  Wu  and  Professor  Marshall  P.  Tulin 


Wave  Generated  At  Bow  of  a  Frigate.  X/Lpp=l.2.  H/X=0.025.  Fn=0.3 


Generated  At  Bow.  A/Lpp«=l.2.  H/X=0.025.  Fn=0.3 


Frigate:  H/X  =  0.025;  Fn  =«  0.30 


elevation,  Ap,  computed  in  time  using  non-linear  slender  body  theory 


tnSuaq  9AQiii  / 


Deck  Wetness,  Slamming,  vs.  Wave  Parameters 


Bow  Wave  Produced  by  a  Frigate  Model 


Bow  Wave  Produced  by  a  Frigate  Model 


Bow  Wave  Produced  by  a  Frigate  Model 


qiSusT]  aAOjft  /  iqSiaH 


l(litne)  /  t(encounter  period) 


|C  Tutiij  ‘ 


Peak  water  elevation,  A^.  computed  In  time  using  non-iinear  slender  body  theory 


■nclosur*  B 
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1  Introduction 

Among  various  problems  of  seakeepiug  qualities  of  ships,  deck  wetness  or  shipping  of  water  becomes 
important  espaciaiiy  in  high  speed  operation  in  rough  seas.  It  is  a  governing  factor  for  determination  of 
freeboard  and  the  design  of  flare  at  the  bow  in  ship  design  practice.  The  existing  method  of  theoretical 
prediction  of  seakeeping  qualities  depends  on  the  linearized  theory  of  the  fluid  motion  around  the  hull. 
In  practice,  hydrodynamic  forces  on  the  hull  are  predicted  by  means  of  the  strip  theory.  It  is  generally 
accepted  that  the  ship  motion  determined  from  the  linear  equation  of  the  motion  with  hydrodynamic 
coefficients,  derived  from  the  liiiearizcrl  theory  in  two  dimensions,  shows  a  fairly  good  prediction  for 
the  ship  beliavior  among  waves,  unless  in  very  severe  conditions.  However  it  is  knosrn  that  some  of 
the  prediction  by  means  of  the  linearized  theory  is  not  necessarily  satisfactory  in  several  problems  of 
seakeepiug  characteristics  including  deck  wetness  and  slamming  impact.  Deck  wetness,  which  is  the 
problem  of  the  present  interest,  is  determined  by  reisttive  height  of  the  deck  above  the  disturbed  sea 
surface.  The  simplest  approximation  by  the  existing  practice  is  that  the  wave  surface  is  taken  as  if  it 
were  undisturbed  by  the  existence  of  the  ship  hull.  It  is  clear  that  this  is  a  too  coarse  assumption,  because 
cIe\'ation  of  the  sea  surface  changes  due  the  disturbance  by  the  hull  to  a  considerable  extent,  so  that  the 
interaction  between  the  hull  and  the  free  surface  should  be  taken  mto  account  by  the  prediction  of  free 
surface  election.  The  heave  of  the  sea  surface  due  to  the  disturbance  by  the  hull  is  called  the  dynamic 
swell.  The  existing  practice  employs  the  strip  theory  with  the  linearized  theory  in  two  dimensions  applied 
to  each  cross  section  of  the  hull.  Takaishi  et  al(1972)  calculated  the  relative  height  of  the  deck  above  the 
disturbed  ware  surface  of  a  contamer  ship  by  the  above  method.  Fairly  good  agreement  was  observed 
between  computation  and  experiment.  Since  this  is  the  case  of  large  container  ship  at  forward  speed 
of  moderate  Froude  number,  the  ship  motion  is  comparatively  mild,  and  the  linearized  theory  seems  to 
work  well.  However  the  condition  is  different  when  the  ship  is  operating  at  high  speed  in  rough  seas. 
Nonlinearit}'  becomes  remarkable  in  the  fluid  motion  around  the  bull,  and  breaking  wave  or  spray  is 
obsen'ed  at  the  bow,  which  is  likely  to  have  a  considerable  effect  on  shipping  water.  Thcie  phenomena 
arc  quite  nonlinear  and  there  has  been  no  reliable  method  so  far  to  predict  deck  wetness  in  such  a  severe 
condition. 

As  the  analytical  solution  is  not  applicable  to  the  nonlinear  boundary  value  problem,  which  appears 
in  the  present  case,  the  analysis  depends  on  the  numerical  method  in  general.  Although  recent  develop¬ 
ment  of  the  computational  fluid  dynamics  enables  the  numerical  solution  of  the  fully  three-dimenaonal 
nonlinear  free  surface  flow  problem,  the  three-dimensional  computation  in  the  present  problem  does  not 
seem  pratical  because  of  large  computer  time  and  insufficient  accuracy.  The  strip  theory  is  the  simplest 
approximation  in  which  nonlinear  effects  can  be  considered.  However  the  strip  theory  is  unable  to  de¬ 
scribe  the  effect  of  forward  speed  to  the  fluid  motion  in  rational  way.  Application  of  the  slender  body 
technique  to  the  ship  moving  through  ambient  ocean  waves  provides  a  rational  approach.  The  existing 
theor}'  of  this  kind  is  based  on  a  scries  expansion  of  the  solution  with  respect  to  two  small  parameters 
which  arc  mutually  independent.  The  slenderness  ratio  of  the  hull  and  the  ratio  of  wave  height  to  the 
ship’s  length  arc  taken  as  the  perturbation  parameters  in  geiicral.  Then  the  flrst  order  solution  is  deter¬ 
mined  from  a  .set  of  linear  boundary  value  problems  in  two  dimensions  in  the  plane  perpendicular  to  the 
longitudinal  axis  of  the  hull.  Thus  the  slender  ship  theory  formulated  in  this  way  belongs  originally  to 
the  liucari'/.cd  approximation.  An  implication  of  the  two  parameter  expansion  is  the  condition  that  the 
aiiiplitiidr  of  the  vertical  inotiu  of  the  hull  as  well  as  the  free  surface  elevation  is  small  as  compared  with 
the  <iimeii.sioii  of  the  crons  .section  of  the  hull,  i.e.  breadth  and  draft.  However  the  motion  amplitude  is 
not  necessarily  small  in  comparison  with  the  draft,  even  though  the  amplitude  of  oscillation  and  wave 
height  arc  l>oth  small  in  conipari.soii  with  the  ship’s  length,  which  is  taken  as  the  reference  length  in  the 
perturbation  .scheme.  This  coiitra«li<;tion  is  the  coii.sequence  of  the  assumption  that  the  slenderness  ratio 
and  the  amplitude  ratio  are  niiitinilly  iudci>eudeiit.  As  a  matter  of  fact,  these  parameters  arc  of  the  same 
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Ollier,  :iii)|  the  se|iiU':iliou  of  two  varitihles  i.s  not  npiilieiiliU*  to  the  present  nntUysis.  Since  tlic  sUnulcrncs.s 
ratio  anil  the  aniplilmle  ratio  are  of  the  .same  order,  it  is  more  ratioaJ  that  they  are  expressed  by  a 
eoiiimon  parameter.  The  perturbation  expansion  with  respect  to  a  single  parameter  Icatb  to  another  set 
of  boundary  value  problems,  which  i.s  iliiferent  from  tlie  existing  .slender  ship  theory  Irased  on  the  two 
parameter  expansion.  A  great  .advantage  of  this  scltcine  is  that  the  first  order  solution  allows  the  large 
ain|>ltlnde  motion  in  the  local  .scale,  even  if  a  small  .amplitude  is  .assumed  in  the  global  sciUc  based  on 
the  ship's  length.  Dilferent  from  the  e.xisting  .slender  ship  theory,  the  boundary  value  problem  becomes 
noiiliear  in  the  new  approach.  This  may  hring  some  complication  in  the  solution  technique.  However  the 
rciiiark;il)lc  jnogress  of  the  coinpntation.al  fluid  dynamics  in  recent  years  h.ss  enabled  the  solution  of  fully 
nonlinear  bonnd.ary  v.alite  itroblem  of  the  free  surface  flotv  to  be  tractable  at  least  iu  two  dimensions. 

.As  far  .as  the  ship  has  no  forward  speed,  the  i)riucip.al  idea  of  the  solution  is  similar  to  the  non¬ 
linear  strip  theory,  by  which  the  fluid  motion  at  each  transverse  plane  is  purely  two-dimensional  and 
independent  each  other.  When  the  finite  forward  speed  is  introduced,  however,  the  fluid  motion  in  a 
transverse  plane  at  one  cross  section  is  subject  to  the  effect  of  the  fluid  motion  at  other  cross  section. 
Then  the  fluid  motion  around  the  hull  becomes  substantially  three-dimensional.  A  great  advantage  of 
the  slender  body  technique  is  that  the  three-dimensional  fluid  motion  is  determined  by  the  solution 
of  the  two-dimensional  boundary  value  problem  at  each  cross  section  even  though  the  fluid  motion  is 
thrcc-dimcnsioual.  Then  the  numerical  solution  for  the  uouiiucar  free  surface  flow  is  much  simplified 
from  the  fully  three-dimensional  solution.  The  numerical  method  for  the  solution  of  fuDy  nonlinear  free 
surface  problem  in  two  dimensions,  which  has  been  developed  recently,  is  available  in  the  solution  of  the 
thrce-dimcusioiial  problem  in  this  manner. 

In  order  to  examine  deck  wetness,  one  must  determine 

1.  Vertical  movement  of  the  deck. 

2.  Height  of  the  free  surface  at  the  hull  surface. 

The  vertical  movement  of  the  deck  is  determined  from  the  equation  of  motion  of  the  hull  Since  the 
fluid  motion  around  the  hull  is  nonlinear,  the  hydrodynamic  forces  and  moments  are  nonlinear  functions 
of  the  ship  and  the  incident  wave.  According  to  experimental  data  however,  the  nonlinearity  of  the 
hydrodynamic  forces  and  moments  acting  on  the  hull  as  a  whole  is  not  remarkable  except  in  such 
extreme  conditions  that  bottom  slamming  takes  place.  It  is  generaUy  accepted  that  the  solution  of 
linear  equations  of  motion  predicts  the  ship  motion  fairly  well,  if  the  hydrodynamic  coeflEidents  in  the 
equatons  are  suitably  chosen.  Correction  may  be  applied  to  the  hydrodynamic  coeflidents  determined 
by  the  linear  strip  theory  in  order  to  illustrate  the  nonlinear  effect.  The  nonlinear  slender  ship  theory 
enables  the  nonlinear  computation  of  hydrodynamic  forces  and  moments.  Then  the  ship  motion  may 
be  determined  by  integration  of  the  nonlinear  equation  of  motion,  but  the  numerical  work  will  become 
enormous.  In  the  present  work,  the  ship  motion  is  assumed  to  be  given  by  a  simple  equation.  Then  the 
boundary  condition  on  the  hull  surface  is  given  by  known  functions,  and  the  fluid  motion  around  the 
hull  is  determined  under  given  boundary  conditions  both  on  the  hull  surface  and  on  the  free  surface. 

Several  problems  are  encountered  in  the  development  of  a  reliable  computation  sdieme.  The  numer¬ 
ical  .solution  shows  instability  sometimes.  Another  problem  of  difficulty  is  treatment  of  the  intersection  of 
the  hull  surface  and  the  free  surface,  when  the  boundary  condition  takes  different  form  on  two  boundaries 
of  different  kinds.  Since  the  level  of  the  point  of  intersection  is  the  measure  of  deck  wetness,  determina¬ 
tion  of  it  is  very  important  in  the  present  problem.  Wave  breaking  and  spray  are  likely  to  appear  on  the 
bow  region.  The  computation  scheme  should  be  able  to  handle  the  fluid  motion  accompanied  by  wave 
breaking  and  spray  generation. 

In  the  present  work,  effort  is  made  to  overcome  the  above  difficulties  in  order  to  find  out  a  method 
of  determining  the  free  surface  elevation  at  the  bow,  which  governs  deck  wetness,  when  wave  breaking 
or  spray  generation  Is  present.  The  aim  of  this  work  is  to  develop  a  method  of  predicting  dedc  wetness, 
which  i.s  more  accurate  tlmii  the  c.xi.stiiig  mothed  based  on  the  linearized  strip  theory. 


2  Outline  of  the  mathematical  formulation 

Take  the  cartesian  coordinates  z,  y,  z  fixed  in  space  with  the  origin  in  the  still  water  plane,  with  z 
and  y-nxes  taken  horizont.olly  and  z-axis  taken  vertically  upwjuds. 

The  fluid  i.s  .xssiiiiied  as  invLscid  and  incompressible.  Then  the  fluid  motion  started  from  rest  is 
irrofatioii.-il  and  i.s  .specilied  by  the  velocity  potential  #  which  .satisfies  the  Laplace  eqntioii  in  the  space 
ocrnpied  by  the  llilid. 


▼ 


|£)  =  0  (1) 

where  sul>scripts  uienu  the  partial  derivativ’cs.  The  boundary  condition  for  ^  consists  of  tn-o  parts.  One 
is  tlic  boundary  condition  on  the  hull  surface  and  the  other  is  on  the  free  surface. 


As  shown  in  Figure  1,  consider  a  ship  moving  through  a  regular  train  of  waves  with  average  forward 
speed  U  in  the  direction  of  positive  x— axis.  Take  another  coordinate  system  (X,K,Z)  fixed  to  the  ship, 
n-itli  .\'-axis  along  the  longitndinal  axis  of  the  ship,  K-axis  athwart  ships  and  Z-axis  in  the  upward 
direction.  The  incident  wave  is  taken  to  propagate  in  the  z-direction  with  the  profile  expressed  by  the 
equation 

z  =  Cy.  s  Asiu(kz  +  («'t  +  So)  (2) 

where  k  s  2-ir/X  s  u*/g,  A  is  the  wave  length  and  ut  is  the  circular  frequency.  The  ship  makes  osdllatious 
with  three  degrees  of  freedom,  i.e.  surge,  heave  and  pitch.  The  motion  of  the  ship  is  expressed  by 
coordinates  of  the  center  of  gravity  (Zo,0,r;„)  in  terms  of  the  z,y,r-axes,  and  the  pitching  angle  $. 
positive  in  the  bow  down  ratation.The  following  transformation  b  valid  between  the  coordinate  system 
(z,y,r)  and  (AM'.Z). 

X  =  (i  -  z„)cos^  -  (z  —  r;,)sm^ 

Y  =  g  (3) 

Z  =  (z-Xo)sin^  +  (z-ro)cosJ 

The  geometry  of  the  hull  surface  b  expressed  by  the  equation 

r  =  F(X,Z)  =  /(x.z,t)  (4) 

The  body  boundary  condition  at  the  hull  surface  Sh  b 

(ifj  ft  +  #*/»  +  #*/»  -  #,  =  0  on  Sh  (a) 

Writing 

Fx  =  dF/dX,  Fz  =  dF/dZ 

one  can  express  (3)  in  the  form 

[H]  (-ig  costf  +  ^  +  (~^o  »>“ ^ 

+#,/,  + =  0  (6) 

where  dots  mean  time  derivatives. 

The  free  surface  condition  consists  of  the  kinematic  and  the  dynamic  conditions.  If  the  free  surface 
b  expressed  by  z  s  (,  the  kinematic  condition  b 

(A'l  Cl  +  -  #.  =  0  atz  =  C  (7) 

The  dynamic  condition  that  the  pressure  b  constant  on  the  free  surface  b 

{D\ 
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(8) 


+  y  { 'f;  4-  +  <l>l)  +  tfC  =  0  at  ;  =  < 

Till*  vi'liirity  |>otru(ial  is  tlic  siiin  of  tlu*  inciiU'iit  wave  |ii>tcutial  0^^.  atul  the  disturbance  potential  <p, 
Stu  ll  as  'P  =  +  0.  We  have  the  eoiulitiun  at  inKiiity  #  =  0„.  at  y/ x’**  +  y'^  -*  oo.  If  the  depth  of 

water  is  inriiiite,  we  have  the  eomlitioii  0.  =  I)  at  r  — •  — oo.  However  these  conditions  for  water  of 
iiifiniti*  streteli  may  not  be  convenient  for  nniiiertcal  solution.  Then  the  width  and  the  depth  of  water 
are  assumed  finite  in  the  actual  computation. 

The  .solution  of  the  boundary  value  problem  defined  above  is  .simplified  to  a  great  extent  by  appli- 
ention  of  the  slender  body  tlieor}*.  The  fiiiulaiiiental  assuiiiption  is  that  the  breadth  D  and  the  draft  d  of 
the  ship  are  both  niiich  smaller  than  the  length  L,  D/L  4^  l,d/L  <ti  1.  Now  let  us  define  the  slenderness 
ratio  as  ;  =  d/L,  and  assume  D/d  =  0(1).  The  slender  body  theory  is  based  on  asingiilar  perturbation 
in  till’  near  field. 

The  formation  is  efficiently  induced  by  the  coordinate  stretching  teclinique.  The  governing  equation 
is  transformed  to  a  system  of  stretched  coordinate  x',  y', such  as  x  =  x',  y  =  ey',  z  =  ez'.  The  solution 
is  e.viiandcd  in  ascending  powers  of  the  governing  equation  is  reduced  to  the  Laplace  equation  in  two 
dimensions  in  the  transverse  y,x— plane, 

m  (9) 

In  applying  .similar  techniques  to  the  boundary  conditions,  we  assume  the  foUowing, 

C/L  =  0{£),  z^/L  =  Ois),  6  =  0{e) 


Then  the  hull  surface  condition  is  transformed  to 

[B\  ^  =  -  [uFx  +  (fffl  +  4  -  Xtf)  Fz]  [1  +  Fl]  (10) 

where  n  is  the  outward  normal  of  the  hull  contour  and  we  can  put  7  =  z  —  Xg.  The  Idnematic  and 
dynamic  free  surface  conditons  become 

[AT  +  =  0  atx  =  C  (11) 

1^1  +  y  +  yC  =  0  at  X  =  f  (12) 

Ill  the  derivation  of  the  above  equations,  we  have  assumed  w\y/djg  =  0{1),  where  ui  is  the  circular 
frequency  of  encounter,  wi  =  w  +  kU,  otherwise  the  free  surface  condition  is  a  trivial. 

Since  the  velocity  potential  is  decomposed  as 


#  =  0„.  +  0 

The  boundary  condition  for  the  disturbance  potential  0  becomes 

[^1  =  -  [uFx  +  (cry  +  4  -  Xy)  Fz]  [i  +  Fl]  ^ 

If  the  incident  wave  is  not  steep,  one  can  write 

3<t>yy  _  0^  ^ 
dz  ~  dt  ^ 

Putting 

v  =  c/y  +  4-Xy-c 

the  hull  surface  condition  is  written  as 

im  ^^=-{UFx  +  VFx){l-\-Fl)-^‘^  on/V 

The  free  surface  elevation  due  to  the  disturbance  is 


(13) 

(14) 

(15) 

(16) 
(17) 


Ci  =  C-Cv  (18) 

Then  the  kinematic  free  surface  condition  can  be  written  as 

(A'l  +  =0  atx  =  C  =  <;r  +  Ci  (19) 

The  dynamic  cumlition  becomes  on  the  other  hand 
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(20) 


PI 


ai  :  ss  C,-  +  Cl 


Ct  iiixl  at  the  free 
initial  condition.  It  is 
above  c([uatioiis  with 
y,  plane  is  defined 


surface  are  determined  by  integration  of  (19)  and  (20)  with  time  with  a  suitable 
more  convenient  to  employ  the  Lagraugean  time  derivative  in  order  to  integrate  the 
time.  The  LagrangeriU  derivative  or  the  material  derivative  in  the  two-dimensional 


by 


£ 

dt 


d  ^  d 


(21) 


Then  the  kinematic  condition  is  written  as 

{A'l  ^  =  on.-  =  C.+(i  (22) 

while  the  dynamic  condition  is  written  as 

PI  j(0j  +  d>I)-«^,Ce-ffCi  at.  =  C.+Ci  (23) 

Thus  the  free  surface  from  =  C  value  of  d>  on  the  &ee  surface  are  determined  from  (22)  and  (23) 
respectively,  but  we  must  know  the  values  of  and  on  the  free  surface  before  integrating  the  above 
equations.  The  fluid  velocities  4>yi4>;  determined  by  the  solution  of  the  two-dimensional  Laplace 
equation 


[L]  ^yy  +  0....  =  0  (24) 

which  satisfies  both  the  hull  surface  condition  (17)  and  the  free  surface  condition  (22)(23). 

It  should  be  noted  that  the  condtion  at  infinity,  ^  0  as  y  -*  d:oo,  is  not  applied,  because  the 

solution  is  v’alid  only  in  the  near  field.  The  solution  of  (24)  is  generally  expressed  in  the  form  -f  ^(z), 
where  4>\  satisfies  the  boundary  conditou  and  vanishes  at  infinity  and  g(x)  is  an  arbitrary  function  of 
X  only,  which  gives  some  indetermiaateness  to  the  solution.  In  order  to  make  the  solution  definite,  one 
has  to  determine  g(x)  by  matching  with  the  far  field  solution  which  may  be  obtained  m  another  way. 
According  to  the  linearized  slender  ship  theory  however,  it  id  proved  that  ^(x}  becomes  negligible  under 
the  condition  u-^d/g  s  0(1).  Therefore  the  condition  at  infinity  holds  in  the  present  case. 


3  Solution  procedure 

The  slender  body  technique  reduces  the  problem  to  the  two-dimensional  boundary  value  problem 
in  the  plane  perpendicular  to  the  longitudinal  axis  of  the  ship  as  explained  in  the  previous  section.  One 
must  remember  that  the  problem  is  formulated  in  the  coordinate  system  fixed  in  space  rather  than  the 
coordinate  system  moving  with  the  ship.  Therefore  the  solution  is  found  in  the  transverse  plane  fixed  in 
space,  and  the  hull  is  moving  relatively  to  this  plane  as  shown  in  Figure  2.  The  cross  section  of  the  hull 
in  the  plane  of  solution  at  one  instant  is  different  from  the  cross  section  in  the  same  plane  at  another 
instant. 
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■rak<'  a  I  l•:lll.sv^•l■sc  ulano  at  a  ('i-rtaiii  iiositiou  r  =  ri  as  tin-  [ilaiii'  of  solution.  The  domain  iu  this 
l<laii*<  o«-i-u|ii<'<l  by  I  ho  lliiid  is  boiiiiilod  by  tUo  hull  cross  soctioii  Fn  and  the  intersection  of  the  free  surface 
Fi- .  The  boiiinlary  value  problem  is  to  .seek  the  solution  of  the  two-diinousional  La|>lacc  equation  (24). 
which  satlslios  the  liull  surfa<  e  condition  (17)  on  Fn  together  with  the  free  surface  condition  (22)  (23) 
on  I'l- .  The  solution  is  obtained  at  each  instant  on  planes  iit  various  positions  x  =  xi.xj.xs,-  •  •.  Then 
the  three-<limcnsiunal  flow  fiehl  is  deternnne<l. 

Tin;  boundary  value  problem  «lescribe<l  above  is  solved  through  two  steps. 

[Step  l| 

A-s-sume  that  the  value  of  on  Fp  is  known  at  a  certain  instant  t  =  ti.  The  normal  velocity  <p„  on 
Fit  i.s  given  by  (17)  if  the  motion  of  the  hull  is  known.  Now  let  us  take  a  closed  domain  circumscribed 
by  T//,  Fp  iuid  a  contour  taken  at  a  great  distance  surrounding  the  fluid  domain,  and  apply  Green’s 
theorem  in  the  domain  fi  inside  this  boundary  contour.  Take  a  field  point  P  inside  fl  and  a  source  point 
Q  on  the  botiutlary  contour  F  =  Fn  +  Fp  +  F^.  Then  Green’s  theorem  gives 

s  '>•  (“) 

where  n  is  the  normal  to  T  in  the  inward  direction  to  fi,  and  subscripts  P,  Q  mean  the  values  at  P  and 
Q  respectively.  If  P  is  a  point  on  F,  the  integral  is  taken  in  the  sense  of  the  Cauchy  principal  value,  and 
the  equation  becomes 

-  -i)’]  (“) 

where  $  =:ir  when  the  contour  has  a  continuous  tangent.  The  normal  velocity  s  d<p/dn  is  given  on 
Fii  by  the  hull  surface  condition  (17)  while  0  is  assumed  to  be  known  on  />.  If  the  contour  ij*,  is  taken 
at  an  infinite  distance,  the  conditon  at  infinity  is  valid  and  the  integral  along  has  no  contribution. 
However  it  is  found  that  taken  at  a  finite  distance  is  more  conveninet  for  numerical  solution.  Here 
we  take  Fx  composed  of  two  vertial  lines  on  both  sides  of  the  hull  and  a  horizontal  line  representing  the 
water  bottom  of  finite  depth,  and  impose  the  conditon  =  0  on  tills  boundary.  The  value  of  <p  oa  Fa 
and  Fx  is  unknown,  while  <(>„  on  Fp  is  unknown.  Then  a  set  of  simidtaneous  integral  equations  for  these 
unknowns  is  derived  from  (26).  The  integral  equation  is  solved  by  a  numerical  method. 

[Step  2] 

The  solution  of  the  integral  equation  obtained  in  Step  1  gives  the  normal  velocity  on  Fp.  Since 
the  velocity  potential  0.is  assumed  to  be  known  on  Fp,  the  tangential  velocity  on  Fp  is  obtained  by 
differentiation  of  ^  along  Fp,  ip,  =  d<t>lds.  The  y  and  z  components  of  the  fluid  velocity  on  the  &ee 
surface  are  determined  by 

=  +  +  (27) 

=  +  (28) 

The  free  surface  elevation  and  the  velocity  potential  on  the  free  surface  at  subsequent  instants  are 
determined  by  the  evolution  equation  (22)  (23)  respectively.  Applying  the  values  of  (27)  and  (28)  to 
(22)  (23),  the  increment  of  Ci  and  in  a  short  interval  At  b  given  by 

40  =  At  dO/dt  (29) 

A<t>  =  At  dtp/dt  (30) 


These  values  arc  calculated  with  respect  to  a  definite  fluid  particle.  The  position  of  the  particle  b 
<lctcrmincd  by 

=  ^  +  Cv  (31) 

Thii.s  the  free  surface  elevation  C  =  Cv  +  Cl  and  the  velocity  potential  on  the  free  surface  at  the  instant 
=  t\+ At  arc  determined,  and  the  boiiiulary  value  problem  in  Step  1  at  the  time  t  =  t]  b  fully  defined. 
The  ste]>  by  .step,  time  mareliiug  procedure  of  Step  1  through  Step  2  determines  the  evolution  of  the 


ilovv  ficUl  in  a  control  t'lanc  at  a  definite  position  x  =  xi.  say.  Since  tlii.s  is  an  initial  value  problem,  the 
initial  condition  at  an  initial  point  must  be  defined.  The  initial  point  is  at  the  intersection  of  the  stem  of 
bow  with  the  free  surface,  and  the  computation  starts  at  the  instant  when  the  initial  point  readies  the 
control  plane.  The  simplest  idea  for  the  initial  condition  is  that  0  =  0  on  the  uudisturlied  free  surface, 
while  the  normal  velocity  on  the  hull  surface  is  given  by  the  hull  surface  condition.  However  it  is  found 
after  some  trial  computation  that  difficulty  appears  in  the  numerical  solution  of  Step  1  by  this  type  of 
initial  condition  because  of  the  singularity  at  the  intersection  of  the  hull  surface  and  the  free  surface. 
The  initial  condition  employed  here  is  as  follows.  In  the  first  time  inter%-al  after  the  initial  point,  the 
scale  of  the  hull  cross  section  in  the  control  plane  is  verj'  small,  that  means  large  Froude  number  of  local 
fluid  motion  and  the  effect  of  gravity  does  not  contribute  to  the  fluid  motion  much.  If  the  ship  has  a 
raked  stem,  the  motion  of  the  hull  contour  in  the  control  plane  is  similar  to  the  vertical  motion  of  a 
sharp  wedge,  for  which  there  is  an  analytical  solution  if  gravity  is  not  present.  Here  we  employ  Mackie’s 
( 1969)  solution  for  the  entry  of  a  sharp  wedge  into  a  free  surface  as  the  initial  conditon.  According  to 
this  theory,  the  velocity  potential  and  the  free  surface  elevation,  when  a  wedge  of  apex  angle  a  enters 
the  free  surface  at  the  velocity  V,  are  given  by 


+  + 


(32) 

(33) 


where  y  and  z  are  normalized  by 


Vt  and  4>  is  normalized  by  V*t. 


Figure  3  Computational  boundary 


4  Numerical  method 

The  integral  equation  in 
is  written  in  the  form 

94>^ 

where 

The  numerical  boundary  is  taken  as  that  in  Figure  3.  On  the  hull  surface  contour  Tm,  is  given  by  the 
boundary  condition,  and  on  tlie  free  surface  contour  />,  ^  is  given  by  the  result  of  computation  in  Step 
2.  On  the  outer  boundary,  wc  set  =  0.  The  outer  boundary  is  divided  into  the  side  wall  Fw  and  the 
w'atcr  bottom  Fg-  The  condition  of  ;nc  horizontal  water  bottom  is  satisfied  by  taking  the  image  with 
respect  to  the  bottom  line  at  z  =  -A.  Hereafter  we  define  Fg  and  Fp  by  including  their  image  contours. 
The  integral  equation  (34)  is  written  for  a  point  P  on  /}/  or  Fw  in  the  form  like 


Step  1  is  solved  by  a  numerical  method.  The  boundary  integral  equation 

0G{P,Q) 


=  /  (p-<^(P>Q)  -  *0 

Jru+rr+r^  \  ^'‘o  ^"o  / 

G{P,  0)  =  y  [(»!.  -  K,)*  +  (-e  -  ip)*] 


(34) 

(35) 
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(36) 


where  C  =  C{P.Q),Cn  =  c>C{P,Q)  /Ohq  for  brevity.  If  the  point  P  is  on  Tf,  the  integral  equation  is 
written  in  the  from 


The  right  hand  sides  of  the  above  equations  are  known  quantities.  The  contour  is  divided  by  discrete 
points  with  small  inter\'als  by  straight  line  segments  connecting  each  adjacent  point.  The  values  of  <t> 
and  are  defined  at  these  points,  and  the  linear  variation  of  ^  and  is  assumed  along  each  segment. 
The  integral  equation  is  then  replaced  by  a  set  of  simultaneous  linear  algebraic  equations.  The  value  of 
0  in  the  integral  equation  in  this  time  is  the  angle  interior  to  the  domain  Cl  between  each  contiguous 
segment. 

The  simultaneous  algebraic  equations  can  be  written  in  matrix  form, 

[[A,]  m]  r  {^}  1  _  \[B,]  (c,i]  f  1 

[iM  m\  \  {#}  /  -  [iB,]  IC,]\  [  {#„}  /  ^38) 

where  mean  values  on  Fh  or  Fw,  and  on  Ff.  If  we  write  the  above  equation  in  simpler  form 

(39) 


the  elements  of  the  matrices  are  given  by 


Aij  —  BiSij  + -  f  aGnds - ^ —  f  Gnds 

~  Jj-i  "*  Jj-i 

rf+x  1  rf+x 

^  ,  i  —  .  I  Gnds - j  sGnds 

^j+X  "■  J)  *i+X  ^  *}  Jj 

Bii  s - i -  r  sGds - /  Gds 

Sj  ~  Sj-l  Jj—l  Sj  —  Sj-l  Jj-l 


...,  rf+i  1  ri+x 

/  Gds - - - /  . 

^i+i  “  Jj  f/+x  ”  Jj 


sGds 


(40) 


(41) 


We  have  a  system  of  linear  algebraic  equations  with  respect  to  tpj  on  Fa  and  Fw  and  on  Ff.  The 
unknown  values  of  <t»  on  Fa  and  Fw  &nd  0  on  />'  are  determined  simply  by  the  matrix  inversion  using 
the  computer.  Thus  the  flow  field  in  a  certain  transverse  plane  T  at  a  certain  instant  (  =  ti  is  completely 


There  are  several  points  which  need  special  treatment.  Since  a  measure  of  deck  wetness  is  relative 
height  of  the  xvare  surface  at  the  hull  surface,  and  the  point  of  intersection  of  the  hull  surface  and  free 
surface  presents  the  top  of  the  wetted  region  on  the  hull  surface,  it  is  obv*ious  that  cautious  treatment 
of  the  intersection  point  is  very  important  for  the  prediction  of  deck  wetness.  Several  researchers  have 
presented  methods  to  treat  this  point,  but  some  of  their  methods  seem  to  be  inaccurate,  and  some  are 
not  suitable  in  the  present  problem.  The  values  of  <p  and  4>n  defined  at  nodal  points  connecting 
adjacent  line  segments  which  replace  the  boundary  contour.  Though  the  potential  is  unique  at  any  nodal 
point,  the  normal  velocity  does  not  have  a  unique  value  at  nodes  where  the  normal  to  the  boundary 
is  not  unique.  In  particular,  the  intersection  point  is  a  special  node  at  which  ^  is  continuous  while 
is  different  on  each  side  because  of  different  boundary  conditions  os  shown  in  Figure  4.  Write  for 
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I  ho  iioriii;il  voUioity  to  tlio  liiill  surfiu-o  atui  0,7  to  tlio  froo  !iiirrai'c  at  the  iiitcracction  point.  Since  0^  is 
known  and  0,7  nukuown,  (39)  can  In*  transforined  to 


(42) 


where  07  €  0„ . . 

Tlie  iinkiiown  07  is  determined  from  tlie  above  ecpiation.  Another  kind  of  nodal  point  which  needs 
.special  treatnient  is  the  tip  of  the  spray  jet  which  .sometimes  occurs.  The  potential  is  given  at  this  point, 
while  the  unknown  normal  velocity  at  the  free  surface  has  different  values  on  both  sides  of  the  jet.  Two 
unkuowii.s  at  the  .same  nodal  point  require  another  crpiation.  The  method  to  solve  this  problem  iin'olves 
discoiitiiiuou.s  elements  iutroducerl  l)y  Drcbbia  and  Douiingttez  (1089)  as  shown  in  Figure  S.  Instead  of 
the  nodal  point  B  connecting  segments  AD  and  DC,  we  take  a  point  B'  on  AD  and  B"  on  DC.  Then 
0  and  0„  arc  calculated  at  B'  and  D".  The  values  of  0„  at  B  on  both  sides  of  the  jet  are  obtained  by 
extrapolation. 

The  next  time  step  is  the  determination  of  the  free  surface  at  the  instant  slightly  later  by  a  short 
time  interval  M.  For  each  time  step,  the  main  computations  described  below  are  performed  for  the 
deformation  of  the  &ee  surface. 

1)  Determination  of  unknowns  0,  #n  :  Equation  (42)  is  rearranged  in  the  form 

“  [iali  -icll]  { {#"} }  *  [irll]  {♦"} 

where  [Px]  ^<i  subvectors  of  [Cj]  in  (42).  The  right  hand  side  of  (43)  is  a  known  vector, 

and  {X}  is  determined  by  the  matrix  inversion. 

2)  Time>stepping  :  The  position  of  fluid  particles  and  values  of  potential  0,  on  the  free  surface  for 
the  next  time  at  t  =  ti  4*  are  determined  by 

y  (ti  +At)  =  y  (ti)  +  ^,(ti)At  (47) 

z  (ti+At)  =  z  (ti)  +  0,(f|)/lt  +  C,(tx+40  -  (^(h)  (48) 

0(ti+4t)  =  0(t,)  +  [i(0*  +  0j)  -  0.Cr  -  (49) 

The  fluid  velocity  0^,  0,  is  calculated  by  (27)  (28),  and  the  tangential  velocity  0,  is  expressed  by 
a  finite  difference  of  0  using  the  three  point  Lagrange  interpolation. 

3)  Smoothing  :  To  depress  a  sawtooth  instability  of  the  wave  profile,  the  five-point  smoothing 
algorithem  is  used  to  filter  the  points  whidi  are  not  equally  spaced  on  the  im  surface,  s,-  in  the 
following  formula  is  the  distance  between  two  nodal  points. 

For  the  first  point  (edge  point)  : 

To  =  -T  -  'M/o  +  /i)  ~  *x(/o  +  A)}  + 

For  the  second  point : 

2(m-  ..) {tTTTI***^*  * ^‘>1 

y-j-^7jj7-[»«(/o  +  /-i)+  •-i(/o  +  /s))} 


(50) 


(51) 


(43) 

(44) 

(45) 

(46) 
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For  the  center  point ; 

/o  =  -rl - “*■  +  ■^i)l  + 

i  ^  s_i  + 


where  i_j  =  5-2  +  5_i,  f  =  ^2  =  +  *2.«‘  =  *i  +  «2  +  *2.”  =  »i  +  <2  +  »3  +«4- 


(32) 


4)  Regridding  :  In  some  cases,  node  points  on  the  free  surface,  especially  near  the  hull,  are  either  too 
close  to  or  too  far  from  each  other,  leading  to  numerical  difiSculties.  The  regridding  technique  is  a 
useful  way  to  keep  the  computation  accurate.  If  the  computed  distance  is  less  than  dmin  or  greater 
than  dmax  which  are  the  input  control  parameters,  the  points  only  near  the  hull  are  redistributed 
by  cubic  spliue  or  linear  interpolation. 


Thus  the  computation  at  the  next  time  step  is  repeated  by  the  above  procedure  1)  through  4). 
Since  the  ship  is  moving  with  forward  speed  U,  the  hull  section  mox’es  astern  in  each  time  inten'al,  and 
readies  the  stern  at  the  time  LjU  after  the  beginning  of  the  above  procedure  at  the  bow.  Therefore  the 
result  of  computation  in  one  round  dose  not  express  the  time  evolution  of  the  free  surface  at  a  definite 
section  of  the  ship  nor  the  free  surface  elevation  in  different  sections  at  the  same  instant.  In  order  to 
complete  the  computation  to  obtain  the  result  for  every  section  throughout  the  length  of  the  ship,  the 
same  computation  procedures  starting  from  the  bow  should  be  repeated  by  shifting  the  initial  instant 
by  the  time  interval  At  =  AxjU,  where  Ax  is  the  interval  of  the  section  for  which  the  result  is  required. 


5  Numerical  Examples 

A  frigate  model  shown  in  Figure  6  is  employed  for  the  numerical  example.  The  computation  of  ship 
motion  in  regular  waves  by  means  of  the  strip  method  and  results  of  towing  tank  experiment  of  thh 
model  have  been  published  by  O’dea  and  Waldeu(1984). 

Numerical  computations  are  carried  out  in  two  conditions.  The  first  is  the  steady  forward  motion 
ill  still  water.  The  Froude  number  varies  from  0.10  to  0.60.  We  take  W  =  60B  and  h  =  1.2L  as  the 
com])utationl  domain.  The  niunbcr  and  intervals  of  the  nodal  points  distributed  on  three  boundaries 
fa,  Tiv  and  ff  arc  as  follows. 

(1)  On  T//  :  The  number  of  nodes  depends  on  the  wetted  depth  of  the  hull  contour,  sudi  as  Ax  =  d/20 
near  the  keel  and  As  s  d/40  near  the  point  of  intersection  with  Ff. 

(2)  On  /tv  :  15  points  are  dustributed  with  an  equal  interval. 
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(;i)  Oil  /'/••  :  150  i»oiiil.s  art’  «listril>iilc<l  from  tin*  iiiti*nnn'tioii  willi  T/;  in  an  cf|iial  interval  As  =  DjlGO, 
anil  Ollier  50  points  are  (li.strilinteil  aloiif;  tin*  roitiaining  jiart  of  the  honinlary  with  arithmetic 
intervals. 

The  niiinorieal  results  at  the  time  interval  At  rIvc  the  free  surface  profile  in  the  transverse  sections 
with  intervals  =  U At.  The  time  step  is  chaseii  in  such  a  way  that  the  niiinerical  results  are  obtained 
at  the  hull  .seetiuns  with  the  intcr\'al  Ax  =  I>/800.  The  reliability  of  the  niimerial  results  is  cotufirmed 
by  ehangiiig  the  intertill  Ax  and  the  grid  dcn.sit.y  on  the  free  surface. 

The  results  of  computation  are  illustrated  in  Figure  7  through  Figure  13.  At  low  Froude  numbers 
less  than  0.2,  generation  of  spray  is  not  observed,  but  diverging  waves  generated  at  the  bow  are  gradually 
steepened  and  eventually  break  at  the  pointed  crest,  suggesting  the  white  cap.  At  Froude  number  0.2, 
overtiiriiiiig  of  the  w.avc  crest  is  observed,  but  the  spray  is  not  developed  yet.  At  Froude  number  0.3, 
generation  of  spray  by  the  overturning  free  surface  is  clearly  observed.  At  Froude  number  0.4,  some 
instability  ai>pcars  oil  the  free  surface  underneath  the  spray  root,  which  may  suggest  the  spilling  type 
of  wave  breaking  instead  of  the  plunging  type.  At  Froude  number  0.5,  some  undulation  appears  on  the 
free  surface  under  the  spray  sheet,  and  this  undulation  grows  up  further  at  Froude  number  0.6.  This 
phenomenon  suggests  another  type  of  wave  breaking. 

The  second  case  Ls  the  same  ship  making  heaving  and  pitching  oscillation  in  sinusoidal  head  sea 
waves  as  expressed  by  eq.(2)  at  uniform  forward  speed  of  Froude  number  0.3.  The  length  of  the  incident 
wave  is  so  chosen  as  X/L  —  1.2,  at  which  the  ma.ximum  amplitude  of  the  oscillations  is  observed  in 
experiment.  The  ship  motion  is  given  by 

Z  =  Z*  sin(a>tt  +  kxg  +  So  +  Ci) 

6  =  sin(a/it  +  kXg  +  Eq  +  Sj) 

where  wi  —  u  +  kU  and  Xg  is  the  x-coordinate  of  the  center  of  gravity  referred  to  the  origin  at  F.P.. 
Values  of  Z^,  Ci  and  £}  are  taken  firom  experimental  data  of  O’dea  et  al. 

The  position  of  the  control  plane  is  defined  by  the  phase  angle  of  the  incident  wave  on  the  control 
plane  at  the  instant  when  the  foremost  point  of  the  bow  arrives  at  this  control  plane  as  shown  in  Figure 
14,  and  the  numerical  procedure  starts  at  this  instant.  If  the  control  plane  fixed  in  space  is  arranged 
with  the  interval  ^xq,  the  phase  difference  between  adjacent  planes  is  Aeo  ~  {ui/U)Axo‘  The  numerical 
results  for  the  free  surface  profile  are  obtained  at  hull  sections  with  the  interval  Ax  =  L/1000,  so  the 
time  iuterv'al  of  the  computation  is  chosen  as  =  0.002256  seconds.  The  reliability  of  the  numerical 
results  is  checked  by  changing  the  interval  Ax  and  the  grid  density  as  before. 

Three  kinds  of  the  wave  steepness  (the  ratio  of  wave  height  to  wave  length)  i.e.  E/X  =  0.02,  0.03 
and  0.04  are  chosen  for  the  numerical  examples.  Computations  are  carried  out  at  control  planes  at 
'0  =  0.0.  ir/4,  ir/2,  3^/4,  x  and  7ir/4.  However  the  preliminary  computation  has  shown  that  there  is  no 
po.ssibility  of  deck  wetness  at  Eq  greater  than  t/2,  so  that  results  with  Eq  =  0.0,  ir/4  and  k/2  will  given 
here. 

1}  Results  otH/X  =0.02 

The  evolution  of  the  free  surface  profile  in  the  bow  region  at  the  control  plane  Eo  =  0.0  is  shown 
in  Figure  15.  At  Section  133{X/L  =  —0.058)  of  the  hull,  the  level  of  the  free  surface  exceeds  the  deck 
height,  but  water  is  pushed  aside  by  the  hull  at  the  deck  side.  The  overturning  of  the  free  surface  results 
the  sheet  of  spray  jet  directed  outward,  so  that  there  is  no  shipping  water  on  board  as  far  as  the  vertical 
movement  of  the  free  surface  relative  to  the  deck  is  upward.  At  Section  201{X/L  =  -0.126)  however, 
the  relative  motion  of  the  free  surface  turns  to  downward.  Then  water  once  lifted  above  the  deck  level 
is  scooped  up  at  the  deck  side,  and  a  small  amount  of  shipping  water  enters  on  dedc  as  shown  in  Figure 
16. 

The  result  at  the  control  plane  Eo  =  x/4  is  shown  in  Figure  17.  The  overturning  free  surface  and  the 
outward  spray  sheet  are  similar  to  the  former  case,  but  shipping  water  resulted  by  the  relative  movement 
of  the  free  surface  as  shown  in  Figure  16  in  the  former  case  does  not  take  place  in  this  case.  Another 
difference  is  the  ap[>earancc  of  undulation  of  the  free  surface  beneath  the  spray  sheet,  which  may  cause 
the  breaking  of  free  surface  behind  this  position. 

The  result  at  the  control  plane  Eq  =  s’/2  is  shown  in  Figure  18.  The  behavior  of  the  free  surface  is 
different  frotii  the  former  cases.  As  the  heave  of  free  surface  cxcecils  the  deck  level  at  Section  S8  much 
clo.ser  to  the  stem  than  in  the  former  cases  due  to  deep  iinmergcncc  of  the  bow,  the  spray  appears  at  the 
tip  of  lio.avcd  water  and  turns  inward  over  the  deck  due  to  outward  movement  of  the  deck  side.  Then 
shipping  water  on  board  is  resulted. 


11 


C(Mii|)iitatiou.s  at  tlio  (.'niitrol  pliuio  fu  =  3ir/4,  ir  aiul  7if/4  ahuw  iio  sliippiiig  water.  As  a  cuuseqiicucc, 
oiir  can  (-oiicIikU'  tliat  dock  wotnoss  begins  to  take  place  in  the  wave  heiglit  H/X  =:  0.02  at  eo  =  O'O 
<t/2.  Tlii.s  ro.sult  is  in  agrcciueiit  with  exporiineiital  data. 

2)  Results  of  H/\  =0.03 

Tlio  behavior  of  the  free  .siirfiice  is  quite  tiilfereut  from  the  former  case,  because  of  the  greater  motion 
of  tlx;  Iiiill.  Heavy  deck  wetness  Ls  reported  in  c.xperimcut.  The  uuiiicrical  result  at  cq  =  0-0  is  shown  in 
Figure  19.  As  the  bow  plunges  into  the  wave  surface,  a  mass  of  water  is  raised  up  on  both  sides  of  the 
deck  in  thick  spray  like  lirpud  walls.  Water  at  the  top  of  the  wall  spreads  in  both  directions,  forming 
.spray  sheet.s  directing  inward  and  outward  to  the  hidl.  The  inward  spray  falls  on  dedc,  resulting  heavy 
deck  wetnes.s.  Tlii.s  plicnomenon  is  similar  to  the  collapse  of  a  liquid  column. 

The  result  at  Jq  =  s-/4  is  shown  in  Figure  20.  .After  the  heaving  free  surface  exceeds  the  deck  side, 
.s|)ray  directing  inward  similar  to  the  case  of  Figure  18  appears  at  the  tip,  bringing  shipping  uniter.  A 
similar  behavior  of  the  free  surface  is  observed  at  the  hull  section  closer  to  the  stem  at  eo  =  ^/-  as 
shown  in  Figure  21. 

According  to  the  experimental  obserxation,  the  bow  of  the  ship  is  lifted  above  the  wave  surface 
almost  clear  of  water  at  the  phase  angle  greater  than  eo  =  a'/2.  Then  slamming  is  likely  to  take  place 
at  the  bow  region.  The  present  computation  scheme  is  not  application  to  this  situation. 

3)  Results  of  H/X  =0.04 

The  result  at  £o  =  0.0  is  shown  in  Figure  22.  The  free  surface  phenomenon  like  the  collapse  of  a 
liquid  column  is  similar  to  that  in  Figure  10,  but  in  a  much  exaggerated  form.  The  spray  sheet  at  the 
top  of  lifted  water  is  stronger,  resulting  heavy  deck  wetness. 

Figure  23  shows  the  result  at  £o  =  ^/A.  The  heaved  free  surface  turns  inboard  as  exceeding  the 
deck  side  and  brings  heavy  shipping  water  in  the  bow  region. 

A  similar  phenomenon  appears  at  the  hull  section  closer  to  the  stem  at  eo  =  */2  as  shown  in  Figure 
24. 


6  Conclusions 

Nonlinear  computations  of  the  free  surface  elevation  at  the  bow  region  of  a  frigate  model  in  the 
steady  forward  motion  as  well  as  in  the  heaving  and  pitching  oscillations  are  carried  out  under  the 
slender  body  approximation. 

The  computation  for  the  steady  forward  motion  well  simulates  the  generation  of  spray  and  the 
breaking  of  bow  waves.  It  is  found  that  the  pattern  of  breaking  wave  changes  according  to  increase  of 
forward  speed. 

The  computation  for  the  oscillating  ship  at  forward  speed  of  Froude  number  0 JO  shows  the  effect 
of  the  incident  wave  height.  At  the  wave  height  of  H/X  =  0.02,  a  mass  of  water  b  raised  above  the  deck 
and  pushed  aside  forming  spray  like  the  plunging  breaker.  However  only  a  small  amount  of  shipping 
water  is  observed.  At  greater  wave  height  H/X  =  0.03  and  0.04,  heavy  deck  wetness  is  resulted  by  falling 
water  from  the  thick  spray  sheet  at  the  lifted  free  surface  over  the  dedc.  Deck  wetness  begins  at  the 
wave  height  H/X  =  0.02,  and  the  result  is  in  good  agreement  with  the  experimental  observation. 
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Figure  7  Bow  wave  breaking  for  a  fixed  ship  in  still  water  at  Fip^O.  10 


Bow  wave  breaking  for  a  fixed  ship  in  still  water  at  Ft$=0.l5 


breaking  for  a  Fixed  ship  in  still  water  at  F/f=0.20 
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Figure  10  Bow  wave  breaking  for  a  fixed  ship  in  still  water  at  F«=0.30 


Figure  1 1  Bow  wave  breaking  for  a  fixed  ship  in  slill  water  at  P'ii=0.40 


Figure  12  Bow  wave  breaking  for  a  fixed  ship  in  still  water  at  Fn=0.50 
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Figure  14  The  position  of  the  No.l  control  plane  (a)  and 
the  expresstion  of  the  initial  phase  of  wave  (b) 
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avc  cicvalion  for  a  .ship  making  heaving  and  pitching  oscillation 
sinusoidal  head  sea  waves  wi(l«  ///A;=().()2  and  e=n/4 
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Figure  18  Wave  elevation  for  a  ship  making  heaving  and  pitching  oscillation 
in  sinusoidal  head  sea  waves  with  and  e=n/2 


'  a  ship  making  heaving  and  pitching  oscillation 
sea  waves  with  ///^0.03  and  c=0.0 
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Wave  elevation  for  a  ship  making  heaving  and  pitching  oscillation 
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hip  making  heaving  and  pitching  oscillation 
waves  with  ///A=0.03  and  t-idl. 
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Figure  24  Wave  elevation  for  a  ship  making  heaving  and  pitching  oscillation 
in  sinusoidal  head  sea  waves  with  and  t-nfl 


